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Abstract 

We present a few-mode Galerkin model for convection in binary fluid layers 
subject to an approximation to realistic horizontal boundary conditions at 
positive separation ratios. The model exhibits convection patterns in form 
of rolls and squares. The stable squares at onset develop into stable rolls at 
higher thermal driving. In between, a regime of a so-called crossroll structure 
is found. The results of our few-mode model are in good agreement with both 
experiments and numerical multi-mode simulations. 
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I. INTRODUCTION 

Compared to convection in ordinary one-component fluids the spatiotemporal properties 
of binary fluids are far more complex. The evolution of the concentration field is governed by 
the interplay of typically strong nonlinear convective transport and mixing, weak dissipative 
solutal diffusion, and the Soret effect [1,2]. The Soret effect generates concentration gradi- 
ents in response to the externally applied temperature difference and to local temperature 
gradients. The strength of the Soret coupling is measured by the dimensionless separa- 
tion ratio ip. The driving mechanisms are therefore controlled by the Rayleigh number R, 
measuring the temperature stress, and by the separation ratio, i.e. the solutal driving. 

In the present paper we focus on 2-dimensional (2D) convective structures consisting 
of straight parallel rolls in one lateral direction, and 3-dimensional (3D) structures that 
look like a nonlinear superposition of two perpendicular roll sets. These structures exist at 
positive separation ratios ip > and arise from a stationary supercritical bifurcation either 
directly out of the ground state or out of a primary convective state. 

At onset, the convection is driven mainly by the solutal gradient established via the Soret 
effect. Therefore this regime is called Soret regime in the literature, see e. g. [3]. The stable 
convection structure is typically a 3D pattern with square symmetry. 

For larger Rayleigh numbers the concentration homogenizes and the fluid behaves more 
like a pure fluid. Convection is now driven mainly by temperature gradients. The stable 
convection pattern is a 2D roll pattern. 

In the intermediate regime, where the stability changes from stationary squares to sta- 
tionary rolls, there exists another 3D structure, the so-called crossroll pattern. This struc- 
ture bifurcates out of the stable square branch and merges with the roll branch at higher 
R. Finally, for slow solutal diffusion, the competition between rolls and squares leads to an 
oscillating behavior in an interval of heating rates around the bifurcation point from squares 
to crossrolls. 

The bifurcation scenario described above has been verifled by several experimental 
groups, e. g. [4,5]. A detailed theoretical insight into the bifurcation scenario has been 
provided by numerical simulations of Ch. Jung et al. [6]. 

The numerical analysis [3,6-8] elucidating the influence of the spatiotemporal behavior 
of the concentration fleld on various properties of convective states for negative and posi- 
tive separation ratio has clearly shown that the success of a model description sensitively 
hinges upon the representation of the concentration fleld. The representation has to capture 
the essence of the spatiotemporal structures following from the combined action of strong 
nonlinear advection and weak diffusion on one hand and the generation of Soret induced 
concentration currents by temperature gradients on the other hand. 

For ip > a. model with few degrees of freedom that reproduces all essentials of the 
bifurcation behavior of the flow amplitude is presently not available. The flrst attempt to 
model this bifurcation topology by Miiller et al. could only generate stable rolls and unstable 
squares [3]. 

Our paper aims at fllling this gap. We present a few-mode Galerkin model which is 
based upon a careful analysis of the concentration balance [3,7-9] in liquid mixtures. With 
it we explain the whole bifurcation scenario from stable squares at onset up to stable rolls. 
The model is an extension of the few-mode model, which was presented in [7] for negative 
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separation ratios. 

We introduce the system and formulate the theoretical task in Sec. II . The main body 
of this paper consists of the two following sections: In Sec. Ill we construct the Galerkin 
model and give a detailed account of how the concentration field is represented. Sec. IV is 
dedicated to a discussion of the results. The convection states are compared in quantitative 
detail with simulations. We summarize our results in Sec. V. 
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II. SYSTEM 

We consider a convection cell of height d. It contains a binary fluid of mean temperature 
T and mean concentration C of the lighter component confined between two perfectly heat 
conducting and impervious plates. This setup is exposed to a vertical gravitational accel- 
eration g and to a vertical temperature gradient AT/d directed from top to bottom. The 
fluid has a density p which varies due to temperature and concentration variations governed 
by the linear thermal and solutal expansion coefficients ct: = ^"^^ ~ ~]>§c' ^^spec- 

tively. Its viscosity is z/, the solutal diffusivity is D, and the thermal diffusivity is n. The 
thermodiffusion coefficient kx quantifies the Soret coupling which describes the driving of 
concentration currents by temperature gradients. 

The vertical thermal diffusion time d"^/ k is used as the time scale of the system and 
velocities are scaled by n/d. Temperatures are reduced by and concentrations by 

,2 

The scale for the pressure is given by Then, the balance equations for mass, momentum, 
heat, and concentration [1,2] read in Oberbeck-Boussinesq approximation [10] 

(2.1) 

aV^u + C7 (T + C) e, (2.2) 

-V ■ [uT - VT] (2.3) 
-V ■ [uC -LV{C- tpT)] . (2.4) 

Here, the currents of heat and concentration, Q and J respectively, are introduced and T 
and C denote deviations of the temperature and concentration fields, respectively, from their 
global mean values T and C. The Dufour effect [10] that provides a coupling of concentration 
gradients into the heat current Q and a change of the thermal diffusivity is discarded in 
(2.3) since it is relevant only in few binary gas mixtures [11]. 

Three parameters enter into the field equations (2.1)-(2.4): the Prandtl number a = i'/k, 
the Lewis number L = D/k, and the separation ratio ^jj = The latter characterizes the 

sign and the strength of the Soret effect. Positive Soret coupling ip induces concentration 
gradients parallel to temperature gradients. In this situation, the buoyancy induced by 
solutal changes in density enhances the thermal buoyancy. When the total buoyancy exceeds 
a threshold convection sets in, typically in the form of squares for positive ip. 

A fourth parameter, the Rayleigh number R = ^^AT measuring the thermal driving 
of the fluid enters the description via the boundary conditions of the temperature fleld (see 
below) . 

The strength of the convection and its influence on convective temperature transport can 
be measured by the Nusselt number = -^(Q ■ ez)xy giving the ratio between the lateral 
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average of the vertical heat current through the system and its conductive contribution. 
In the basic state of quiescent heat conduction its value is thus 1 and larger than 1 in all 
convective states. 

Solving the partial differential equations (2.1)-(2.4) requires boundary conditions for the 
fields. We use realistic no slip conditions for the top and bottom plates at z = 

1 , 

u{x,y,z = ±^;t) = , 
and assume perfectly heat conducting plates by 

T{x,y,z = ±^]t) = T^R . 
Furthermore, impermeability for the concentration is guaranteed by 

e, ■ J = -Ld, (C - ijT) {x,y,z = ±^;t) = . (2.5) 



III. MODE SELECTION AND GALERKIN MODEL 
A. Temperature and velocity field 

To describe the convective state we use the Galerkin method. We should stress that 
we restrict ourselves to the description of extended patterns that are periodic in the lateral 
directions x and y with a certain lateral periodicity length A = 27r/k and fixed phases. We 
take k = 3.117, i.e. the critical wave number of the pure fluid. 

The temperature field which consists of a linear conductive profile —z and a convective 
deviation is truncated as 

T{x,y, z;t) = —Rz+Too2{t)V2sm{27rz) 

+ [Tioi(t) cos(A;x) + Toii(t) cos{ky)] Vq.cos^tiz) (3.1) 
+Tii2{t) cos{kx) cos{ky)V2sm{27rz) . 

The indices (/, m, n) of the amplitudes denote the expansion in x-, y- and ^-direction, 
respectively. If we restricted ourselves to 2D convective patterns only, e. g. homogeneous in 
x-direction, only the modes with the first index equal to zero would be taken into account, 
reducing the temperature ansatz to that of a Lorenz model. In order to have an analogous 
representation of roll patterns in y-direction we also include the mode Tqh. As we see in 
Eq. (3.1) a further nonlinear mode, T112, caused by the interaction of the two roll patterns 
is also introduced. The contribution of this mode vanishes for 2D patterns. Note that only 
modes with an even index sum I + m + n appear in the expansion, which is due to a mirror 
glide symmetry of the patterns studied [9]. 

Choosing the first Chandrasekhar function, Ci{z) [8,12] as an approximation to the z- 
profile of the critical mode, the vertical component w of the velocity field reads 



w{x,y,z,t) = [wioi{t) cos{kx) + woii{t) cos{ky)]Ci{z) 



(3.2) 
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The X- and ?/-components can be derived by using the incompressibihty criterion (2.1). As 
the velocity field is dominated by its critical modes not only at onset but also far beyond, 
these two modes suffice to describe the patterns to be studied. When discussing bifurcation 
diagrams we will use the amplitudes wioi and wqu as order parameters. 



B. Approximation of the boundary condition 

To select adequate concentration modes a detailed analysis of the concentration balance 
and of the field structure of roll and square states is necessary. 
By introducing the combined field ( = C — ipT expanded as 

oo oo oo 

C(x, I/, 2, ^) = Cimn cos{lkx) cos(m/c?/)sc„(2) , (3.3) 

1=0 m=0 n=0 

it is possible to fulfill the realistic boundary condition for the concentration current using 
appropriate trigonometric functions for the sc„(2;) to ensure d^C = at the plates. 

But as we aim to construct a few-mode model, we have to consider the balance equations 
using the concentration field instead of the C-field. As discussed in [7,9,13] for stationary 
and travelling patterns in detail, the reason lies in the balance equation for the ^-field: 

(^t + V ■ u) C = - tpV^T . (3.4) 

The fiuid mixtures which we refer to have Lewis numbers up to 2 orders of magnitude smaller 
than the separation ratio. In this work we consider L = 0.01, tp = 0.1 and L = 0.0045, 
= 0.23 for example. Thus, in order to approximate Eq. (3.4) consistently we have to take 
into account temperature modes which are up to 2 orders of magnitude smaller than the (- 
modes. Thus, despite the fact that the temperature field is well described by the critical and 
the first nonlinear mode alone, many more modes would be required when using Eq. (3.4). 
To avoid this problem, we expand the concentration field itself. As a consequence, we can 
guarantee the boundary condition (2.5) only approximatively. 



1. Approach 

We adopt Hollinger's successful approximation for the boundary condition, which he pro- 
posed in [9]. His investigations on the concentration current induced by advection, diffusion 
and the Soret coupling support to demand impermeability only for the lateral mean of the 
concentration field. In the lateral mean Eq. (2.5) becomes 

{e,-3),y = -Ld,Ciatiz)~LRtpN = . (3.5) 

Here 

Ciat{z) = (C(x, y, z; t))^y (3.6) 

denotes the lateral mean of the concentration field and is the Nusselt number introduced 
in Sec. II. This suggests the following multi-mode Galerkin ansatz for the concentration 
field 
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C{x,y, z;t) = —RipNz + y > > cimnCos{lkx) cos{mky)sCn{z) . 




(3.7) 



;=0 m=0 n=0 



In order to avoid introducing temperature modes in Eq. (3.7) via we approximate the 
boundary condition (3.5) by dz{C{x, y, z; t))xy{H/'2) = —Rip. This deviates from the correct 
value by a factor equal to the Nusselt number = 0(1). This approximation can be 
understood as the leading term in an amplitude expansion of around the conductive state 
with N = 1. 



In this subsection we present a brief comparison of the results of an "exact" multi-mode 
Galerkin simulation [8] which fulfills the boundary conditions of the concentration field 
exactly by introducing the C-held given by Eq. (3.3), with another multi-mode simulation. 
For the latter, we use the same field ansatz for temperature and velocity as for the exact 
one but the concentration field ansatz is given by Eq. (3.7), with replaced by 1. We refer 
to this simulation as using approximated boundary conditions (ABCs) as opposed to exact 
boundary conditions (EBCs). 

Fig. la shows the behavior of the vertical mean Cvcrt(3;) of the concentration for roll 
solutions orientated in |/-direction for a reduced Rayleigh number r = R/Rc{ip = 0) = 
i?/1708 = 1.5 where the homogenization of the concentration outside of narrow boundary 
layers is already apparent. Both, the EBC and the ABC simulation generate practically the 
same result; the ABCs have vanishing influence on the vertical mean. 

The vertical dependence of the first lateral Fourier mode 



displayed in Fig. lb shows also good agreement in the bulk. The influence of the different 
boundary conditions is restricted to the boundary layer only. The ABC solution has a 
vanishing slope at the plates as has any contribution to the concentration field except for 
the vertical mean. The EBC solution ends with a finite slope at smaller values. 

These results imply that the ABCs, that allow the construction of few-mode models, 
suffices to describe the concentration field in the bulk. 



Analogous to the description of the temperature and the velocity field we want to in- 
clude only few modes describing the concentration field. Still, we found that 10 modes are 
necessary to reproduce the bifurcation scenario. The concentration field ansatz which we 
finally used for the few-mode model is given by 



2. Exact versus approximate concentration boundary condition 




A 



(3.8) 



C. Selecting the concentration field modes 



C{x,y,z;t) = —Ripz + Cooi{t)V2sm{7iz) + Coo3{t)\/2sm{37iz) 
+ [ciooit) cos{kx) + coio(t) cos(%)] 
+ [cio2(^) cos(A;x) + coi2{t) cos{ky)] cos(27r2;) 



(3.9) 
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+ [c2oi(^) cos(2/i;x) + Co2i{t) cos{2ky)] -\/2 sin(7r2) 
+ [c3oo(^) cos(3A;x) + co3o(t) cos(3A;y)] . 

For all amplitudes the index sum / + m + n is odd, reflecting the mirror glide symmetry 
already mentioned above. To check and justify the selection of these modes we analyze 
the lateral and vertical mean of the concentration field in the following subsection and thus 
demonstrate the relevance of these modes. We compare our results with the multi-mode 
simulations presented in [8]. Since we found similar results for both, square and roll solutions, 
we will concentrate on rolls here and do not present the detailed comparison of the square 
results. 

1. Lateral average of the concentration field 

The lateral mean Ciat(-z) of the concentration field is determined by modes of the form 
coon- Our few-mode ansatz takes two such modes into account, namely cqoi and cqos- Fig- 
ure 2a shows Ciat(^) for two different reduced Rayleigh numbers r. For r = 1 (thin solid 
line), i.e., at the boundary between Soret and Rayleigh region, the homogenization of Ciat(^) 
in the bulk with pronounced concentration boundary layers at the lower and upper plate 
becomes already apparent. The profile hardly changes when switching to a higher r = 1.5 
(thin dashed line). The two modes, cqoi and cqos, that are also taken into account describe 
this behavior well (thick lines), whereas one Fourier mode alone would be insufficient. 

2. Critical modes of the concentration field 

The critical modes have a lateral dependence cos (^^) depending on the orientation of 
the rolls studied and are thus linear combinations of modes with amplitude f*^""). For each 
case, we again took two modes into account, ('^^''°) and (^^°^). The vertical variation of the 

' ^ ' Vcoio/ Vcoi2/ 

first lateral Fourier mode Ciq{z) for rolls orientated in y-direction is shown in Fig. 2b. For 
both, r = 1 and r = 1.5, the field in the bulk is reproduced well, whereas the boundary 
layers are not properly resolved. As we saw before, the reason for this does not lie in an 
insufficient number of modes but already in the ABC as the ABC multi-mode simulation 
has shown the same deficiency. 

3. Vertical average of the concentration 

Cvertix), the vertical mean of C{x,y, z;t) is represented by the first two lateral modes 
of the roll patterns, sin and sin (3^^) . In Fig. 2c the vertical average for the multi- 
mode simulation is compared to the results of our few-mode model. The multi-mode result 
shows concentration peaks between the convection rolls and homogenization inside them. 
The inclusion of the higher harmonics in the few-mode model allows to reflect this to some 
extent, at least for r = 1 where the concentration peaks are still relatively broad. The more 
narrow peaks at r = 1.5 are not resolved anymore. However, even then there is still good 
agreement further away from the peaks. 
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4- Further modes 

The modes with lateral variation sin (g^^) are decoupled from the others without the 

inclusion of further nonlinear modes. Adding the modes c^2oi^ oc cos (2^^) sin(7r2;) allows for 

driving them via these new modes and W(un\ in the nonlinearity of the concentration field 

equation (2.4). Studies of smaller models indicate that stable squares at onset are impossible 
without including these modes [3,13]. 

Similarly, the modes C(oi2\ discussed in Sec. IIIC2 serve to drive the mode C003. While 

U02J 

this mode would still be driven by W(ioi\ and cnoo\ due to the different 2;-expansion of w 

\oii) \oio) 

and C, the coupling is weak compared to the one via C|'oi2\. 

{102) 

IV. RESULTS 

In this section, we elucidate the roll, square, and crossroll solutions of our model. First, 
we discuss three different bifurcation scenarios. Then, in Sec. IV B, we focus on the oscillat- 
ing patterns generated in the few-mode model in contrast to the oscillating crossroUs which 
appear in the exact simulations. Finally, in Sec. IV C, we present the phase diagram of the 
few-mode model and discuss it in the light of numerical results. 

A. Bifurcation scenario 

In Fig. 3, results for the few-mode model are plotted for three different Lewis numbers in 
the interesting range of heating rates where the transition between small-amplitude stable 
squares in the Soret regime and higher-amplitude stable rolls in the Rayleigh regime takes 
place. All plots presented were calculated at ^ = 0.23 and a = 27. These parameters can 
be realized by ethanol-water mixtures and have also been studied in [6]. 

The upper plot shows a bifurcation diagram representing the three different stationary 
structures found. The square (S) branch is denoted by squares, the roll (R) branch by circles. 
Each is represented by a single curve: For square structures it is w^qi = uigii while for rolls 
one of the amplitudes is zero, say, and thus only Wiqi is plotted. It is it'ioi/'iffoi ~ 
throughout the diagram, or in other words, the sum 1^101 +'W^oii approximately the same for 
both structures. For the third structure, the stationary crossrolls (CR), the two amplitudes 
w^f^ and are different and nonzero, and are thus represented by two curves marked by 
up and down triangles respectively. The crossrolls bifurcate out of the square branch with 
equal amplitude but their difference grows with growing r until wff^ meets the roll branch 
and becomes zero. The plot has been calculated for L = 0.003 but it does not change 
qualitatively and hardly quantitatively for the values of L considered below. 

The second plot shows the real eigenvalues tied to the stationary instabilities for the 
same value of L. They are again identified by square, circle, and triangle symbols. On the 
small r side rolls exhibit a positive eigenvalue and are thus unstable. The square eigenvalue 
is negative even though this is hardly visible on the scale displayed; the smallness of the 
eigenvalue is a consequence of the very slow concentration dynamics (L -C 1) in the Soret 
regime. On the large r side the signs of the eigenvalues are reversed: rolls are now stable 
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while squares are not. At the bifurcation points where the square eigenvalue crosses the zero 
axis the crossrolls appear and vanish again where the roll eigenvalue crosses the zero axis. In 
between, crossrolls are the only stationary stable structure and transfer this stability from 
the squares to the rolls. 

The stationary bifurcation properties of the few-mode Galerkin model, and in particular 
the sequence of stable squares, stable crossrolls, and stable rolls agree well with experimental 
and theoretical results for the system [4-6] for not too small L. However, at very small L an- 
other, time-dependent crossroU pattern appears in the vicinity of the bifurcation point from 
squares to stationary crossrolls. The time-dependent crossrolls emerge from an oscillatory 
perturbation and are thus represented by a complex eigenvalue. 

The real parts of the most important complex eigenvalues for L = 0.003 and two further 
Lewis numbers in our few-mode model are presented in the lower three plots of Fig. 3. The 
results show that when L is small enough, an oscillatory perturbation destroys the stability 
of the square pattern already before the bifurcation point where the stationary crossrolls 
appear. These crossrolls are then also oscillatory unstable and gain stability only at higher 
r before meeting the roll branch (cases L = 0.003 and L = 0.0037). The real part of the 
crossroU eigenvalue has a local minimum at r ^ 0.96 such that crossrolls might temporarily 
gain stability before becoming unstable again. This happens for L = 0.0037. They might 
also lose stability to an oscillatory perturbation only later on while being stable at the 
bifurcation point from the squares. This is the case for L = 0.0045. For even larger L the 
stationary crossrolls remain stable against oscillatory perturbations everywhere (not shown). 
Rolls can become unstable against oscillatory perturbations too, but we only found this to 
be the case at heating rates below the crossroll-roll bifurcation point where they are already 
stationary unstable. 

B. Oscillations 

In [6], Ch. Jung et al. describe the behavior of oscillating crossrolls, as they appear in 
their multi-mode simulations at small L between the regimes of stable squares and stable 
stationary crossrolls. In this structure, the leading amplitudes u^ioi and wou oscillate around 
one of the square fixed points with opposite phase. With growing r the oscillation becomes 
increasingly anharmonic until the oscillatory state disappears in a subharmonic bifurcation 
cascade, that is associated with an entrainment process. 

We studied time-dependent patterns in our few-mode model in the r-range where all 
three stationary structures are unstable. We used the same time-integration method as in 
[6]. While the oscillatory instabilities of our few-mode model occur in the correct region 
of parameter space, the oscillating patterns of the model differ from the oscillating cross- 
rolls discussed in [6]. We found different oscillatory regimes like those presented in Fig. 4 
calculated for a Lewis number of 0.0037. 

The pattern at r = 0.95 is a 2D pattern in which only Wioi depends on time, while 
Won remains at zero (or vice versa). For r = 0.98, u;ioi and wou are equal, preserving the 
square symmetry. This pattern oscillates around one of the square fixed points with growing 
amplitude, followed by chaotic flips between the vicinities of the two square fixed points. 

Since these time-dependent structures disagree qualitatively with those found in the 
multi-mode simulations, we did not study them further. We conclude that a few-mode 
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model does not suffice to capture the nature of the oscillatory crossroUs. 

C. Phase diagram 

Figure 5 illustrates the phase diagrams of the few-mode model in 2D planes of the 4D 
parameter space of r, L, and a. In the following, the results in Fig. 5 are compared to 
the EBC results (Fig. 6) taken from [6]. 

We consider first the L — r plane in Fig. 5a. Qualitatively, the phase diagram is very 
similar to the full numerical EBC result. For small enough L the conductive state loses its 
stability against stable squares which in turn lose stability at higher r against oscillating 
structures discussed in the last subsection. The oscillations are then replaced by stationary 
crossrolls. Those finally merge with the roll branch. For L > 0.005 the oscillatory regime is 
absent, and squares lose their stability directly to stationary crossrolls. At L ^ 1, outside 
of the plotted interval rolls are stable at onset. Note the dent in the region of oscillations at 
L f» 0.004, r 0.96 where the pattern sequence becomes more complicated. Three example 
cases have been studied in Sec. IV A. 

Fig. 6a looks qualitatively very similar. Only the dent in the oscillatory crossroll region 
is absent and the sequence of patterns with increasing r is for small L always squares - 
oscillatory crossrolls - stationary crossrolls - rolls. However, there are two main quantitative 
differences. First, the region of oscillations is larger in the multi-mode simulation; the direct 
squares - stationary crossrolls transition happens only at L > 0.02. The second important 
quantitative difference between multi-mode and few-mode model lies in the position of 
the phase space interface between roll and stationary crossroll patterns. The maximum in 
Fig. 5a lies at r = 1.27, whereas in Fig. 6a the maximum of this curve lies at r = 1.79 . In [6], 
the authors point out that the data gained in numerical simulations become quantitatively 
unreliable for small L combined with large r. Comparing finite difference calculations and 
Galerkin models with different spatial resolutions they found that the range of existence of 
crossrolls shrinks with decreasing resolution which also explains the location of the crossroll- 
roU boundary in our few-mode model. 

Next, we want to study the ^/'-dependence in Fig. 5 and Fig. 6 . Note that our chosen 
fixed value of L = 0.0045 is the one where the dent in the region of oscillations is located in 
Fig. 5a and where the above described sequence of patterns does not occur with growing r. 
Thus, the results in Fig. 5b look qualitatively different for this L value from those of Fig. 6b. 
Stable stationary crossrolls exist above the oscillatory structures as well as below for all ijj 
where oscillatory structures exist at all. This was already seen for the special case = 0.23 
which is again marked by a dotted line. 

Finally, let us consider the a — r plane in Fig. 5c. The stationary stability curves are 
independent of a. The reason lies in the ansatz for the velocity. Our ansatz (3.2) takes only 
the critical modes into account, which eliminates the nonlinearity in Eq. (2.1) and leads to 
cr-independent stationary solutions. In the multi-mode diagram in Fig. 6, only the region of 
stable squares remains roughly independent of a whereas the range of existence of oscillatory 
and stationary crossrolls shrinks with decreasing a. Again, there are qualitative differences 
to Fig. 6c concerning the appearance of oscillations. Oscillatory regimes are found in two 
different regions of Fig. 5c and do not exist at small a. 
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To summarize our comparison between few-mode model and multi-mode simulations, we 
found similar phase diagrams in the L — r plane. However, there are qualitative differences 
for the ip — r and a — r plane concerning the locations of oscillatory crossroUs due to the 
fact that the Lewis number L = 0.0045 studied in [6] lies within a small range of L-values 
where the pattern sequence is more complicated in the few-mode model. 

Finally, we would like to point at another comparison between multi-mode simulations 
using ABCs and EBCs in Fig. 6b: We plotted the results of the ABC multi-mode model 
in the ip — r plane (dashed lines). Compared to the EBC results the ranges of existence 
for all occurring patterns are nearly the same. Only the phase space boundary between 
rolls and crossroUs that sensitively hinges on the model size is shifted significantly. There is 
no qualitative difference in the location of the oscillatory region as found in the few-mode 
model. These results underline again the usefulness of the ABCs especially considering the 
fact that stability properties are in general more model-sensitive than fixed point properties. 

V. CONCLUSION 

This paper presents results of a new few-mode Galerkin model for Rayleigh-Benard 
convection in binary mixtures with positive separation ratio. With only 16 modes it is 
able to produce bifurcation and phase diagrams similar to those found in experiment or 
stemming from much more extensive simulations. In particular, it reproduces a sequence 
of stable squares, crossroUs, and roUs. The success of this model is based on and is due to 
a carefully chosen concentration field ansatz. This ansatz was obtained from comparisons 
with multi-mode simulations. 

We used an approximation to the impermeability boundary condition allowing us to keep 
the number of temperature field modes at a minimum. The concentration field ansatz then 
fulfills the impermeability condition at the plates in a strict sense only for the quiescent 
conductive state. Comparisons between results of multi-mode Galerkin simulations for rolls 
and squares using exact and approximated boundary conditions show that both lead to 
essentially the same concentration field in the bulk. Significant differences appear only in 
the concentration boundary layers. Further comparisons show that the ABC treatment 
also leads to qualitatively the same phase diagram. This is an interesting result as former 
theoretical models with permeable boundary conditions lead to very different features. It 
appears that the impermeability condition for the conductive state is the most important 
condition to fulfill. 

In the few-mode model, we found stable squares and unstable rolls bifurcating out of 
the quiescent conductive state for a wide range of parameters. A transition to stable rolls 
at higher temperature differences (Rayleigh regime) happens via a crossroU pattern. This 
pattern bifurcates forward out of the squares and merges with the rolls at larger r. Previous 
few-mode models also show an intermediate regime between Soret and Rayleigh regime, but 
the behavior of the system was different from the simulations. 

The results for the phase diagrams prove the success of our model, as we have found 
qualitatively similar stability domains for the stationary structures. While our model fails 
to generate oscillating crossroUs as they appear in multi-mode simulations [6], similar oscil- 
latory structures have been found in roughly the same parameter region. 
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Our model is analytically manageable, and further studies can possibly pave the way for 
a better understanding of the behavior and the role of the physically important modes. 
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FIG. 1. Comparison of multi-mode EBC simulations to multi-mode ABC simulations of roll 
convection, (a) Vertical mean Cvcrti^^) of the concentration field versus x. (b) First lateral Fourier 
mode Ciq{z) versus z. Parameters are ip = 0.1, a = 10, L = 0.01, and r = 1.5. 
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FIG. 2. The concentration field of the few-mode model (thick lines) compared to EBC numer- 
ical calculations (thin lines) at r = 1 (solid) and r = 1.5 (dashed) for rolls, (a) Lateral mean, (b) 
Lateral Fourier mode Cio{z). (c) Vertical mean. Parameters are tp = 0.1, cr = 10, L = 0.01. 
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FIG. 3. Bifurcation diagrams and corresponding eigenvalues of our 16 modes model, (a) bifur- 
cation diagram of the velocity modes, (b) important real eigenvalues for (a), (c-e) real parts of the 
important complex eigenvalues. Square (roll) branches are symbolized by squares (circles). Thick 
lines and triangles mark the crossroll structure. Parameters are a = 27, ip = 0.23. 
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FIG. 4. Oscillations in the few-mode model. The dashed lines mark the fixed points of the 
stationary square patterns. Parameters are a = 27, ip = 0.23, and L = 0.0037. 
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FIG. 5. Phase diagrams of stable convective states in the few-mode model. The dotted lines 
refer to the parameters L = 0.0045, ifj = 0.23 and a = 27, respectively. Two of them are kept fixed 
in each of the diagrams while the third is varied. The ordinate scale is logarithmic for r < 0.8 to 
show the bifurcation threshold rstat{k = 3.117; V'j-Z^) of the conductive state. 
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